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ABSTRACT 

We use optical data on 10 Kuiper Belt objects (KBOs) to investigate their rotational properties. Of 
the 10, three (30%) exhibit light variations with amplitude Am > 0.15 mag, and 1 out of 10 (10%) 
has Am > 0.40 mag, which is in good agreement with previous surveys. These data, in combination 
with the existing database, are used to discuss the rotational periods, shapes, and densities of Kuiper 
Belt objects. We find that, in the sampled size range, Kuiper Belt objects have a higher fraction 
of low amplitude lightcurves and rotate slower than main belt asteroids. The data also show that 
the rotational properties and the shapes of KBOs depend on size. If we split the database of KBO 
rotational properties into two size ranges with diameter larger and smaller than 400 km, we find that: 
(1) the mean lightcurve amplitudes of the two groups are different with 98.5% confidence, (2) the 
corresponding power-law shape distributions seem to be different, although the existing data are too 
sparse to render this difference significant, and (3) the two groups occupy different regions on a spin 
period vs. lightcurve amplitude diagram. These differences are interpreted in the context of KBO 
collisional evolution. 

Subject headings: Kuiper Belt objects — minor planets, asteroids — solar system: general 



1. INTRODUCTION 

The Kuiper Belt (KB) is an assembly of mostly small 
icy objects, orbiting the Sun beyond Neptune. Kuiper 
Belt objects (KBOs) are lik ely to be remnants of outer 
solar system planetesimals l|Jewitt fe Luulfl993fl . Their 
physical, chemical, and dynamical properties should 
therefore provide valuable information regarding both 
the environment and the physical processes responsible 
for planet formation. 

At the time of writing, roughly 1000 KBOs are known, 
half of which have been followed for more than one op- 
position. A total of « 10 5 objects larger t han 50 km are 
thoug ht to orbit the Sun beyond Neptune l|Jewitt fc Luul 
2000) . Studies of KB orbits have revealed an intricate dy- 
namical structure, with signatures of interactions with 
Neptune (Malhotra 1995). The size distribution fol- 
lows a dif ferential power-law o f index q — 4 for bodies 
> 50km l|Truiillo et al.M2001aD, becoming s lightly shal- 
lower at smaller sizes ((Bernstein et al.ll2004)) . 

KBO colou rs show a large diversity, from slightly blue 
to very red llLuu fc Jewittl Il99fil iTegler fc Romanishinl 
120001 Ijewitt fe Luull200H) . and seem to cor relate with in- 
clinat i on and/or perihelion di s tance ( e.g.. Uewitt fc Lu~ 
I 2001L iDoressoundiram et all 120021 iTruiillo fc B rown 
2002). The few low- resolution optical and near-IR KBO 
spectra are mostly featureless, with the exception of a 
weak 2 iv ca. water ice absorption line present in some 
of them HBrown et a,UlT999l Uewitt fc Luul 12001. and 
strong methane absorption on 2003 UB313 ijBrown et al.l 
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HI). 

About 4% of known KBOs are binarie s with separa- 
tions larger than 0'.'15 ijNoll et alJ 12002^ . All the ob- 
served binaries have primary-to-secondary mass ratios 
ps 1. Two bina r y crea tion models have been proposed. 
IWeidenschillind l|2002|) favours the idea that binaries 
form in three-body encounters. This model requires a 
100 times denser Kuiper Belt at the epoch of binary for- 
mation, and predicts a higher abundan ce of large sepa- 
ratio n binaries. An alternative scenario l)Goldreich et alJ 
2002), in which the energy needed to bind the orbits of 
two approaching bodies is drawn from the surrounding 
swarm of smaller objects, also requires a much higher 
density of KBOs than the present, bu t it predicts a larger 
fractio n of close binaries. Recently, Shc ppard fc JewitH 
( 2004) have shown evidence that 2001 QG 2 g 8 could be a 
close or contact binary KBO, and estimated the fraction 
of similar objects in the Belt to be ~ 10%-20%. 

Other physical properties of KBOs, such as their 
shapes, densities, and albedos, are still poorly con- 
strained. This is mainly because KBOs are extremely 
faint, with mean ap parent red magnitude mji ^23 
ijTruiillo et al.ll2001blL 

The study of KBO rotational properties through time- 
series broadband optical photometry has proved to be the 
most successful technique to date to investigate some of 
these physical properties. Light variations of KBOs are 
believed to be caused mainly by their aspherical shape: 
as KBOs rotate in space, their projected cross-sections 
change, resulting in periodic brightness variations. 

One of the best examples to date of a KBO lightcurve 
- and what c an be learned from it - is that of 
f20000') Varuna l|Jewitt fc Sheppardll2002T) . The authors 
explain the lightcurve of (20000) Varuna as a consequence 
of its elongated shape (axes ratio, a/b ~ 1.5). They fur- 
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Fig. 1. — Frame-to-frame photometric variances (in magnitudes) 
of all stars (gray circles and black crosses) in the (35671) 1998 SN165 
(a) and (38628) Huya (b) fields, plotted against their relative mag- 
nitude. The trend of increasing photometric variability with in- 
creasing magnitude is clear. The intrinsically variable stars clearly 
do not follow this trend, and are located towards the upper 
left region of the plot. The KBOs are shown as black squares. 
(35671) 1998 SN155, in the top panel shows a much larger variabil- 
ity than the comparison stars (shown as crosses, see Section 13.11 . 
while (38628) Huya is well within the expected variance range, 
given its magnitude. 

ther argue that the object is centripctally deformed by 
rotation because of its low density, "rubble pile" struc- 
ture. The term "rubble pile" is generally used to refer 
to gravitationally bound aggregates of smaller fragments. 
The existence of rubble piles is thought to be due to con- 
tinuing mutual collisions throughout the age of the solar 
system, which gradually fracture the interiors of objects. 
Rotating rubble piles can adjust their shapes to balance 
centripetal acceleration and self-gravity. The resulting 
equilibrium shapes have been studied in the extreme case 
of fluid b odies, and depend on the body's density and 
s pin rate lChandrasekharlll969|) . 

lLacerda fc Luul l|2003l hereafter LL03a) showed that 
under reasonable assumptions the fraction of KBOs 
with detectable lightcurves can be used to constrain 
the shape distribution of these objects . A follow-up 
ijLuu fc Lacerdal 1200.11 hereafter lLL 03b) on this work, 
using a database of lightcurye properties of 33 KBOs 
(jSheppard fc JewittJ l200l 12003). shows that although 
most Kuiper Belt objects (~ 85%) have shapes that are 
close to spherical (a/b < 1.5) there is a significant frac- 
tion (~ 12%) with highly aspherical shapes (a/b > 1.7). 

In this paper we use optical data on 10 KBOs to inves- 
tigate the amplitudes and periods of their lightcurves. 
These data are used in combination with the existing 
database to investigate the distributions of KBO spin 
periods and shapes. We discuss their implications for 
the inner structure and collisional evolution of objects in 
the Kuiper Belt. 

2. OBSERVATIONS AND PHOTOMETRY 

We collected time-series optical data on 10 KBOs at 
the Isaac Newton 2.5m (INT) and William Herschel 4m 
(WHT) telescopes. The INT Wide Field Camera (WFC) 
is a mosaic of 4 EEV 2048x4096 CCDs, each with a 



pixel scale of O'.'33/pixel and spanning approximately 
Ilf3x22f5 in the plane of the sky. The targets are im- 
aged through a Johnson R filter. The WHT prime focus 
camera consists of 2 EEV 2048x4096 CCDs with a pixel 
scale of O'.'24/pixel, and covers a sky-projected area of 
2x8.'2xl6.'4. With this camera we used a Harris R fil- 
ter. The seeing for the whole set of observations ranged 
from 1.0 to 1.9"FWHM. We tracked both telescopes at 
sidereal rate and kept integration times for each object 
sufficiently short to avoid errors in the photometry due 
to trailing effects (see Table QJ. No light travel time cor- 
rections have been made. 

We reduced the data using standard techniques. The 
sky background in the flat-fielded images shows varia- 
tions of less than 1% across the chip. Background vari- 
ations between consecutive nights were less than 5% for 
most of the data. Cosmic rays were re moved with the 
package LA-Cosmic l|van Dokkumll2CJ0l . 

We performed aperture photometry on all objects 
in the field using the SExtractor software package 
ijBertin fc Arnoutsl 11996(1 . This software performs cir- 
cular aperture measurements on each object in a frame, 
and puts out a catalog of both the magnitudes and the 
associated errors. Below we describe how we obtained a 
better estimate of the errors. We used apertures ranging 
from 1.5 to 2.0 times the FWHM for each frame and se- 
lected the aperture that maximized signal-to-noise. An 
extra aperture of 5 FWHMs was used to look for pos- 
sible seeing dependent trends in our photometry. The 
catalogs were matched by selecting only the sources that 
arc present in all frames. The slow movement of KBOs 
from night to night allows us to successfully match a 
large number of sources in consecutive nights. We dis- 
carded all saturated sources as well as those identified to 
be galaxies. 

The KBO lightcurves were obtained from differential 
photometry with respect to the brightest non-variable 
field stars. An average of the magnitudes of the brightest 
stars (the "reference" stars) provides a reference for dif- 
ferential photometry in each frame. This method allows 
for small amplitude brightness variations to be detected 
even under non-photometric conditions. 

The uncertainty in the relative photometry was calcu- 
lated from the scatter in the photometry of field stars 
that are similar to the KBOs in brightness (the "com- 
parison" stars, see Fig^l. This error estimate is more 
robust than the errors provided by SExtractor (see be- 
low), and was used to verify the accuracy of the latter. 
This procedure resulted in consistent time series bright- 
ness data for ~ 100 objects (KBO + field stars) in a time 
span of 2-3 consecutive nights. 

We observed Landolt standard stars whenever condi- 
tions were photometric, and used them to calibrate the 
zero point of the magnitude scale. The extinction coeffi- 
cient was obtained from the reference stars. 

Since not all nights were photometric the lightcurves 
are presented as variations with respect to the mean 
brightness. These yield the correct amplitudes and peri- 
ods of the lightcurves but do not provide their absolute 
magnitudes. 

The orbital parameters and other properties of the 
observed KBOs are given in Table M Tables El H 
and list the absolute i?-magnitude photometric mea- 
surements obtained for (19308) 1996 T0 66 , 1996 TS 66 , 
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Fig. 2. — Stacked histograms of the frame-to-frame variance (in magnitudes) in the optical data on the "reference" stars (in white), 
"comparison" stars (in gray), and the KBO (in black). In c), e), and j) the KBO shows significantly more variability than the comparison 
stars, whereas in all other cases it falls well within the range of photometric uncertainties of the stars of similar brightness. 

We begin by investigating if the observed brightness 
variations are intrinsic to the KBO, i.e., if the KBO's 
intrinsic brightness variations are detectable given our 
uncertainties. This was done by comparing the frame- 
to-frame scatter in the KBO optical data with that of 
(~ 10 — 20) comparison stars. 

To visually compare the scatter in the magnitudes of 
the reference stars (see Section[2Jl, comparison stars, and 
KBOs, we plot a histogram of their framc-to-frame vari- 
ances (see Fig. |2J). In general such a histogram should 
show the reference stars clustered at the lowest variances, 



(35671) 1998 SNi 65 , and (19521) Chaos, respectively. Ta- 
bles and [S] list the mean-subtracted i?-band data for 
(79983) 1999 DF 9 and 2001 CZ 31 . 

3. LIGHTCURVE ANALYSIS 

The results in this paper depend solely on the ampli- 
tude and period of the KBO lightcurves. It is therefore 
important to accurately determine these parameters and 
the associated uncertainties. 

3.1. Can we detect the KBO brightness variation? 
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TABLE 1 

Observing Conditions and Geometry 



Object Designation 


ObsDate a 


Tel. b 


Seeing c 


MvtRt d 


ITime e 


RA £ 




dec g 




R h 


A 1 


a 1 










f" 

L 




["/hr] 


[sec] 


[hhmmss] 


[<""] 




[AU] 


[AU] 


[deg] 


(19308) 


1996 T0 66 


99-Oct-01 


WHT 


1 


8 


2.89 


500 


zo 


•j.j 


Iv) 


tuo 


90 


42 


45.950 


44.958 


0.1594 




1996 TS 6 6 


99-Sep-30 


WHT 


1 


3 


2.62 


400,600 


uz 


9f : ; 


uu 


191 


II 


03 


38.778 


37.957 


0.8619 




1996 TS 66 


99-Oct-01 


WHT 


1 


1 


2.67 


600 


02 


26 


02 


lOI 


40 


50 


38.778 


37.948 


0.8436 




1996 TS 66 


99-Oct-02 


WHT 


1 


5 


2.70 


600,900 


02 


25 


58 


4-91 


40 


35 


38.778 


37.939 


0.8225 


(35671) 


1998SNi65 


99-Sep-29 


WHT 


1 


5 


3.24 


360,400 


23 


32 


46 


—01 


18 


15 


38.202 


37.226 


0.3341 


(35671) 


1998SNi 65 


99-Sep-30 


WHT 


1 


4 


3.22 


360 


23 


32 


41 


—01 


18 


47 


38.202 


37.230 


0.3594 


(19521) 


Chaos 


99-Oct-01 


WHT 


1 





1.75 


360,400,600 


03 


44 


37 


191 


30 


58 


42.399 


41.766 


1.0616 


(19521) 


Chaos 


99-Oct-02 


WHT 


1 


5 


1.79 


400,600 


m 

u-> 


41 


91 


191 

-|-Z1 


9U 


54 


42.399 


41.755 


1.0484 


(79983) 


1999 DFg 


Ol-Feb-13 


WHT 


1 


7 


3.19 


900 


10 


27 


04 


+09 


15 


16 


39.782 


38.818 


0.3124 


(79983) 


1999 DF 9 


Ol-Feb-14 


WHT 


1 


6 


3.21 


900 






OU 


1 no 
Tuy 


1 

1U 


25 


39.783 


38.808 


0.2436 


(79983) 


1999 DFg 


Ol-Feb-15 


WHT 


1 


4 


3.22 


900 


10 


26 


46 


+09 


46 


50 


39.783 


38.806 


0.2183 


(80806) 


2000 CM105 


01-Feb-ll 


WHT 


1 


5 


3.14 


600,900 


09 


18 


48 


+19 


41 


59 


41.753 


40.774 


0.1687 


(80806) 


2000 CM105 


Ol-Feb-13 


WHT 


1 


4 


3.12 


900 


09 


18 


39 


+19 


42 


40 


41.753 


40.778 


0.2084 


(80806) 


2000 CM105 


Ol-Feb-14 


WHT 


1 


5 


3.11 


900 


09 


18 


34 


+19 


43 


02 


41.753 


40.781 


0.2303 


(66652) 


1999 RZ 253 


01-Sep-ll 


INT 


1 


9 


2.86 


600 


22 


02 


57 


-12 


31 


06 


40.963 


40.021 


0.4959 


(66652) 


1999 RZ253 


Ol-Sep-12 


INT 


1 


1 


2.84 


600 


22 


02 


53 


-12 


31 


26 


40.963 


40.026 


0.5156 


(66652) 


1999 RZ253 


Ol-Sep-13 


INT 


1 


8 


2.82 


600 


22 


02 


49 


-12 


31 


49 


40.963 


40.033 


0.5381 


(47171) 


1999 TC 3 e 


01-Sep-ll 


INT 


1 


9 


3.85 


600 


00 


16 


49 


-07 


34 


59 


31.416 


30.440 


0.4605 


(47171) 


1999 TC 36 


Ol-Sep-12 


INT 


1 


1 


3.86 


900 


00 


16 


45 


-07 


35 


33 


31.416 


30.437 


0.4359 


(47171) 


1999 TC36 


Ol-Sep-13 


INT 


1 


8 


3.88 


900 


00 


16 


39 


-07 


36 


13 


31.416 


30.434 


0.4122 


(38628) 


Huya 


Ol-Feb-28 


INT 


1 


5 


2.91 


600 


13 


31 


13 


-00 


39 


01 


29.769 


29.021 


1.2725 


(38628) 


Huya 


01-Mar-01 


INT 


1 


8 


2.97 


360 


13 


31 


09 


-00 


38 


23 


29.768 


29.009 


1.2479 


(38628) 


Huya 


Ol-Mar-03 


INT 


1 


5 


3.08 


360 


13 


31 


01 


-00 


36 


59 


29.767 


28.987 


1.1976 




2001 CZ31 


01-Mar-01 


INT 


1 


3 


2.72 


600,900 


09 


00 


03 


+16 


29 


23 


41.394 


40.522 


0.6525 




2001 CZ31 


Ol-Mar-03 


INT 


1 


4 


2.65 


600,900 


08 


59 


51 


+16 


30 


01 


41.394 


40.539 


0.6954 



a UT date of observation 
^Telescope used for observations 
c Avcragc seeing of the data [ ] 
^Average rate of motion of KBO ["/hr] 
integration times used 

*Right ascension 
^Declination 
^KBO— Sun distance 

^BO-Earth distance 

■'Phase angle (Sun-Object-Earth angle) of observation 



TABLE 2 
Properties of Observed KBOs 



Object Designation 


Class a 


H b 


i c 




a c 






[mag] 


[deg] 




[AU] 


(19308) 1996 TOee 


C 


4.5 


27.50 


0.12 


43.20 


1996 TS 66 


C 


6.4 


7.30 


0.13 


44.00 


(35671) 1998SNi 65 


c f 


5.8 


4.60 


0.05 


37.80 


(19521) Chaos 


C 


1.9 


12.00 


0.11 


45.90 


(79983) 1999 DF 9 


C 


6.1 


9.80 


0.15 


46.80 


(80806) 2OOOCM105 


C 


6.2 


3.80 


0.07 


42.50 


(66652) 1999RZ253 


C 


5.9 


0.60 


0.09 


43.60 


(47171) 1999 TC 36 


Fb 


1.9 


8.40 


0.22 


39.30 


(38628) Huya 


P 


1.7 


15.50 


0.28 


39.50 


2001 CZ31 


C 


5.4 


10.20 


0.12 


45.60 



a Dynamical class (C — Classical KBO, P — Plutino, b — binary 
KBO) 

Absolute magnitude 
c Orbital inclination 
^Orbital eccentricity 
°Scmi-major axis 

^This object as a classical-type inclination and eccentricity but 
its semi-major axis is much smaller than for other classical KBOs 



followed by the comparison stars spread over larger vari- 
ances. If the KBO appears isolated at much higher vari- 
ances than both groups of stars (e.g., Fig. then its 
brightness modulations are significant. Conversely, if the 



TABLE 3 
Photometric measurements of 
(19308) 1996 T0 66 









UT Date a 


Julian Date b 


[mag] 



1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 
1999 Oct 



1.84831 
1.85590 
1.86352 
1.87201 
1.87867 
1.88532 
1.89302 
1.90034 
1.90730 
1.91470 



2451453. 
2451453. 
2451453. 
2451453. 
2451453. 
2451453. 
2451453. 
2451453. 
2451453. 
2451453. 



34831 
35590 
36352 
37201 
37867 
38532 
39302 
40034 
40730 
41470 



21.24 
21.30 
21.20 
21.22 
21.21 
21.28 
21.27 
21.30 
21.28 
21.31 



0.07 
0.07 
0.07 
0.07 
0.07 
0.07 
0.06 
0.06 
0.06 
0.06 



a UT date at the start of the exposure 
Julian date at the start of the exposure 

c Apparcnt red magnitude; errors include uncer- 
tainties in relative and absolute photometry added 
quadratically 



KBO is clustered with the stars (e.g. Fig- I2f ) , any peri- 
odic brightness variations would be below the detection 
threshold. 

Figure shows the dependence of the uncertainties on 
magnitude. Objects that do not fall on the rising curve 
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TABLE 4 

Photometric measurements of 1996 TSee 



UT Date a 


Julian Date b 




[mag] 




1999 


Sep 


30. 


.06087 


O A P" 1 A P" "1 

2451451 


.56087 


21 


.9 1 


± 





.U i 


1999 


Sep 


30 


1 \ i • 1 ' ill! 

.06628 


2451451 


.56628 


21 


.83 


± 





,U / 


1999 


Sep 


30. 


07979 


i a tr i a pt i 

2451451 


.57979 


21 


.76 


± 





IY7 

u * 


1999 


Sep 


30 


.08529 


o a tr 1 a pt i 

2451451 


.58529 


A 1 


.71 


± 


0. 


,U / 


1999 


Sep 


30 


.09068 


o a tr 1 a pt i 

2451451 


.59068 


21 


.75 


± 


0. 


U / 


1999 


Sep 


30 


.09695 


O A P* "1 A P" "1 

2451451 


.59695 


21 


.67 


± 





U / 


1999 


Sep 


30. 


.01250 


O A P" 1 A P" "1 

2451451 


.60250 


O "1 

21 


.77 


± 





IY7 

u * 


1999 


Sep 


30 


.10936 


2451451 


.60936 


Z 1 


.76 


± 


0. 


,UD 


1999 


Sep 


30. 


.11705 


2451451 


.61705 


-il 


.80 


± 





,uo 


1999 


Sep 


30. 


.12486 


o j p" i /i p" I 

2451451 


.62486 


21 


.77 


± 





,UD 


1999 


Sep 


30 


.13798 


n /I CT1 /I PT 1 

2451451 


a i too 

.63798 


■yi 

Z 1 


.82 


± 


0. 


,U / 


1999 


Sep 


30 


.14722 


2451451 


.64722 


■yi 
21 


.74 


± 





rift 
,UD 


1999 


Sep 


30. 


.15524 


o J p" i /i p" I 

2451451 


.65524 


O "1 

21 


.72 


± 





,UD 


1999 


Sep 


30 


.16834 


n /I CT1 /I PT 1 

2451451 


.66834 


■yi 
A 1 


.72 


± 


0. 


,Uo 


1999 


Sep 


30 


.17680 


2451451 


.67680 


21 


.83 


± 





U i 


1999 


Sep 


30. 


.18548 


2451451 


.68548 


O "1 
21 


.80 


± 





,UD 


1999 


Sep 


30. 


.19429 


2451451 


.69429 


21 


.74 


± 





IY7 

u * 


1999 


Sep 


30 


.20212 


O /I CT1 /I PT 1 

2451451 


.70212 


21 


.78 


± 


0. 


,U / 


1999 


Sep 


30 


.21010 


2451451 


.71010 


■yi 
21 


.72 


± 





,U / 


1999 


Sep 


30. 


.21806 


O J P" 1 /I P" "1 

2451451 


.71806 


-il 


.76 


± 





no 

.uy 


1999 


Sep 


30. 


.23528 


2451451 


.73528 


A 1 


.73 


± 





,U / 


1999 


Sep 


30. 


.24355 


2451451 


.74355 


O "1 

21 


.74 


± 





Uo 


1999 


Oct 


01 


.02002 


2451452 


.52002 


O "1 

21 


.81 


± 





,UD 


1999 


Oct 


01 


.02799 


2451452 


.52799 


21 


.82 


± 





.06 


1999 


Oct 


01 


.03648 


2451452 


.53648 


21 


.81 


± 





.06 


1999 


Oct 


01 


.04422 


2451452 


.54422 


21 


.80 


± 





.06 


1999 


Oct 


01 


.93113 


2451453 


.43113 


21 


.71 


± 





.06 


1999 


Oct 


01 


.94168 


2451453 


.44168 


21 


.68 


± 





.06 


1999 


Oct 


01 


.95331 


2451453 


.45331 


21 


.73 


± 





.06 


1999 


Oct 


01 


.97903 


2451453 


.47903 


21 


.69 


± 





.06 


1999 


Oct 


01 


.99177 


2451453 


.49177 


21 


.74 


± 





.06 


1999 


Oct 


02 


.00393 


2451453 


.50393 


21 


.73 


± 





.05 


1999 


Oct 


02 


.01588 


2451453 


.51588 


21 


.78 


± 





.05 


1999 


Oct 


02 


.02734 


2451453 


.52734 


21 


.71 


± 





.05 



TABLE 5 
Photometric measurements of 
(35671) 1998SNi 65 



m R 



a UT date at the start of the exposure 
k Julian date at the start of the exposure 
c Apparent red magnitude; errors include uncertainties 
in relative and absolute photometry added quadratically 



traced out by the stars must have intrinsic brightness 
variations. By calculating the mean and spread of the 
variance for the comparison stars (shown as crosses) we 
can calculate our photometric uncertainties and thus de- 
termine whether the KBO brightness variations are sig- 
nificant (>3<t). 

3.2. Period determination 

In the cases where significant brightness variations 
(see Section I3.1|l were detected in the lightcurves, the 
phase dispersion m inimization method was used (PDM, 
Stcllingwerl ll978|) to look for periodicities in the data. 
For each test period, PDM computes a statistical param- 
eter 6 that compares the spread of data points in phase 
bins with the overall spread of the data. The period that 
best fits the data is the one that minimizes 8. The advan- 
tages of PDM are that it is non-parametric, i.e., it does 
not assume a shape for the lightcurve, and it can handle 
unevenly spaced data. Each data set was tested for pe- 
riods ranging from 2 to 24 hours, in steps of 0.01 hr. To 
assess the uniqueness of the PDM solution, we generated 
100 Monte Carlo realizations of each lightcurve, keeping 
the original data times and randomizing the magnitudes 
within the error bars. We ran PDM on each generated 
dataset to obtain a distribution of best-fit periods. 



UT Date a 


Julian Date b 




[mag] 




1999 


Sep 


29. 


.87384 


2451451 


.37384 


21. 


.20 


± 


0. 


06 


1999 


Sep 


29. 


.88050 


2451451 


.38050 


21. 


.19 


± 


0. 


05 


1999 


Sep 


29. 


.88845 


2451451 


.38845 


21 


.18 


± 


0. 


05 


1999 


Sep 


29. 


.89811 


2451451 


.39811 


21 


.17 


± 


0. 


05 


1999 


Sep 


29. 


.90496 


2451451 


.40496 


21 


.21 


± 


0. 


05 


1999 


Sep 


29. 


.91060 


2451451 


.41060 


21 


.21 


± 


0. 


05 


1999 


Sep 


29. 


.91608 


2451451 


.41608 


21. 


.18 


± 


0. 


05 


1999 


Sep 


29. 


.92439 


2451451 


.42439 


21. 


.25 


± 


0. 


05 


1999 


Sep 


29. 


.93055 


2451451 


.43055 


21. 


.24 


± 


0. 


05 


1999 


Sep 


29. 


.93712 


2451451 


.43712 


21. 


.26 


± 


0. 


06 


1999 


Sep 


29. 


.94283 


2451451 


.44283 


21. 


.25 


± 


0. 


06 


1999 


Sep 


29. 


.94821 


2451451 


.44821 


21 


.28 


± 


0. 


06 


1999 


Sep 


29. 


.96009 


2451451 


.46009 


21 


.25 


± 


0. 


06 


1999 


Sep 


29. 


.96640 


2451451 


.46640 


21 


.21 


± 


0. 


06 


1999 


Sep 


29. 


.97313 


2451451 


.47313 


21 


.17 


± 


0. 


05 


1999 


Sep 


29. 


.97850 


2451451 


.47850 


21 


.11 


± 


0. 


05 


1999 


Sep 


29. 


.98373 


2451451 


.48373 


21 


.12 


± 


0. 


06 


1999 


Sep 


29. 


.98897 


2451451 


.48897 


21. 


.15 


± 


0. 


06 


1999 


Sep 


29. 


.99469 


2451451 


.49469 


21 


.15 


± 


0. 


06 


1999 


Sep 


29. 


.99997 


2451451 


.49997 


21. 


.16 


± 


0. 


.06 


1999 


Sep 


30 


.00521 


2451451 


.50521 


21 


.12 


± 


0. 


.06 


1999 


Sep 


30. 


.01144 


2451451 


.51144 


21 


.09 


± 


0. 


.06 


1999 


Sep 


30 


.02164 


2451451 


.52164 


21 


.18 


± 


0. 


06 


1999 


Sep 


30. 


.02692 


2451451 


.52692 


21 


.17 


± 


0. 


.06 


1999 


Sep 


30. 


.84539 


2451452 


.34539 


21. 


.32 


± 


0. 


.06 


1999 


Sep 


30 


.85033 


2451452 


.35033 


21. 


.30 


± 


0. 


06 


1999 


Sep 


30 


.85531 


2451452 


.35531 


21. 


.28 


± 


0. 


.06 


1999 


Sep 


30. 


.86029 


2451452 


.36029 


21. 


.31 


± 


0. 


.06 


1999 


Sep 


30. 


.86550 


2451452 


.36550 


21. 


.21 


± 


0. 


.06 


1999 


Sep 


30 


.87098 


2451452 


.37098 


21 


.26 


± 


0. 


06 


1999 


Sep 


30. 


.87627 


2451452 


.37627 


21. 


.28 


± 


0. 


06 


1999 


Sep 


30. 


.89202 


2451452 


.39202 


21. 


.23 


± 


0. 


05 


1999 


Sep 


30 


.89698 


2451452 


.39698 


21. 


.30 


± 


0. 


.06 


1999 


Sep 


30 


.90608 


2451452 


.40608 


21. 


.20 


± 


0. 


05 


1999 


Sep 


30. 


.91191 


2451452 


.41191 


21. 


.26 


± 


0. 


.05 


1999 


Sep 


30. 


.92029 


2451452 


.42029 


21. 


.15 


± 


0. 


05 


1999 


Sep 


30. 


.92601 


2451452 


.42601 


21. 


.19 


± 


0. 


05 


1999 


Sep 


30. 


.93110 


2451452 


.43110 


21 


.11 


± 


0. 


05 


1999 


Sep 


30. 


.93627 


2451452 


.43627 


21. 


.16 


± 


0. 


.05 


1999 


Sep 


30. 


.94858 


2451452 


.44858 


21. 


.18 


± 


0. 


05 


1999 


Sep 


30. 


.95363 


2451452 


.45363 


21. 


.16 


± 


0. 


05 


1999 


Sep 


30. 


.95852 


2451452 


.45852 


21. 


.13 


± 


0. 


05 


1999 


Sep 


30. 


.96347 


2451452 


.46347 


21 


.17 


± 


0. 


.05 


1999 


Sep 


30. 


.96850 


2451452 


.46850 


21 


.16 


± 


0. 


05 


1999 


Sep 


30. 


.97422 


2451452 


.47422 


21 


.18 


± 


0. 


05 


1999 


Sep 


30. 


.98431 


2451452 


.48431 


21 


.18 


± 


0. 


05 


1999 


Sep 


30. 


.98923 


2451452 


.48923 


21 


.17 


± 


0. 


.06 


1999 


Sep 


30. 


.99444 


2451452 


.49444 


21 


.16 


± 


0. 


05 


1999 


Sep 


30. 


.99934 


2451452 


.49934 


21 


.20 


± 


0. 


05 


1999 


Oct 


01 


.00424 


2451452 


.50424 


21 


.16 


± 


0. 


.05 


1999 


Oct 


01 


.00992 


2451452 


.50992 


21 


.18 


± 


0. 


.06 



a UT date at the start of the exposure 
b Julian date at the start of the exposure 
c Apparent red magnitude; errors include uncertainties 
in relative and absolute photometry added quadratically 



3.3. Amplitude determination 

We used a Monte Carlo experiment to determine the 
amplitude of the lightcurves for which a period was 
found. We generated several artificial data sets by ran- 
domizing each point within the error bar. Each artificial 
data set was fitted with a Fourier series, using the best-fit 
period, and the mode and central 68% of the distribution 
of amplitudes were taken as the lightcurve amplitude and 
la uncertainty, respectively. 

For the null lightcurves, i.e. those where no significant 



6 
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TABLE 6 

Photometric measurements of (19521) Chaos 



UT Date a 


Julian Date b 




[mag] 




1999 


Oct 


01 


.06329 


2451452 


.56329 


20 


.82 


± 





.06 


1999 


Oct 


01 


.06831 


2451452 


.56831 


20 


.80 


± 





.06 


1999 


Oct 


01 


.07324 


2451452 


.57324 


20 


.80 


± 





.06 


1999 


Oct 


01 


.07817 


2451452 


.57817 


20 


.81 


± 





.06 


1999 


Oct 


01 


.08311 


2451452 


.58311 


20 


.80 


± 





.06 


1999 


Oct 


01 


.08801 


2451452 


.58801 


20 


.76 


± 





.06 


1999 


Oct 


01 


.09370 


2451452 


.59370 


20 


.77 


± 





.06 


1999 


Oct 


01 


.14333 


2451452 


.64333 


20 


.71 


± 





.06 


1999 


Oct 


01 


.15073 


2451452 


.65073 


20 


.68 


± 





.06 


1999 


Oct 


01 


.15755 


2451452 


.65755 


20 


.70 


± 





.06 


1999 


Oct 


01 


.16543 


2451452 


.66543 


20 


.72 


± 





.06 


1999 


Oct 


01 


.17316 


2451452 


.67316 


20 


.72 


± 





.06 


1999 


Oct 


01 


.18112 


2451452 


.68112 


20 


.71 


± 





.06 


1999 


Oct 


01 


.18882 


2451452 


.68882 


20 


.73 


± 





.06 


1999 


Oct 


01 


.19652 


2451452 


.69652 


20 


.70 


± 





.06 


1999 


Oct 


01 


.20436 


2451452 


.70436 


20 


.69 


± 





.06 


1999 


Oct 


01 


.21326 


2451452 


.71326 


20 


.72 


± 





.06 


1999 


Oct 


01 


.21865 


2451452 


.71865 


20 


.72 


± 





.06 


1999 


Oct 


01 


.22402 


2451452 


.72402 


20 


.71 


± 





.06 


1999 


Oct 


01 


.22938 


2451452 


.72938 


20 


.72 


± 





.06 


1999 


Oct 


01 


.23478 


2451452 


.73478 


20 


.71 


± 





.07 


1999 


Oct 


01 


.24022 


2451452 


.74022 


20 


.72 


± 





.07 


1999 


Oct 


02 


.04310 


2451453 


.54310 


20 


.68 


± 





.06 


1999 


Oct 


02 


.04942 


2451453 


.54942 


20 


.69 


± 





.06 


1999 


Oct 


02 


.07568 


2451453 


.57568 


20 


.74 


± 





.07 


1999 


Oct 


02 


.08266 


2451453 


.58266 


20 


.73 


± 





.06 


1999 


Oct 


02 


.09188 


2451453 


.59188 


20 


.74 


± 





.06 


1999 


Oct 


02 


.10484 


2451453 


.60484 


20 


.75 


± 





.06 


1999 


Oct 


02 


.11386 


2451453 


.61386 


20 


.77 


± 





.06 


1999 


Oct 


02 


.12215 


2451453 


.62215 


20 


.77 


± 





.06 


1999 


Oct 


02 


.13063 


2451453 


.63063 


20 


.78 


± 





.06 


1999 


Oct 


02 


.13982 


2451453 


.63982 


20 


.79 


± 





.06 


1999 


Oct 


02 


.14929 


2451453 


.64929 


20 


.71 


± 





.07 



TABLE 7 

Relative photometry measurements of 
(79983) 1999 DF 9 



a UT date at the start of the exposure 
b Julian date at the start of the exposure 
c Apparent red magnitude; errors include uncertainties 
in relative and absolute photometry added quadratically 



variation was detected, we subtracted the typical error 
bar size from the total amplitude of the data to obtain 
an upper limit to the amplitude of the KBO photometric 
variation. 

4. RESULTS 

In this section we present the results of the lightcurve 
analysis for each of the observed KBOs. We found sig- 
nificant brightness variations (Am > 0.15 mag) for 3 out 
of 10 KBOs (30%), and Am > 0.40 mag for 1 out of 10 
(10%). This is consistent with p reviously publis hed re- 
sults: iShermard Jewitti (J2p02, hereafter ISJOl found 
a fraction of 31% with Am > 0.15 mag and 23% with 
Am > 0.40 mag, both consistent with our results. The 
other 7 KBOs do not show detectable variations. The 
results are summarized in Table 

4.1. 1998 SN m5 

The brightness of (35671) 1998 SNi 65 varies signifi- 
cantly (> 5a) more than that of the comparison stars 
(see Figs. H and Et)- The periodogram for this KBO 
shows a very broad minimum around P = 9hr (Fig. [3^). 
The degeneracy implied by the broad minimum would 
only be resolved with additional data. A slight weaker 
minimum is seen at P = 6.5 hr, which is close to a 24 hr 
alias of P = 9 hr. 



UT Date a 


Julian Date b 




5m r c 
[mag] 




2001 


Feb 


13. 


.13417 


2451953 


.63417 


0. 


21 


± 





.02 


2001 


Feb 


13. 


.14536 


2451953 


.64536 


0. 


20 


± 





.03 


2001 


Feb 


13. 


.15720 


2451953 


.65720 


0. 


17 


± 





.0 1 


2001 


Feb 


13. 


.16842 


2451953 


.66842 


0. 


06 


± 





.03 


2001 


Feb 


13. 


.17967 


2451953 


.67967 


-0. 


08 


± 





.02 


2001 


Feb 


13. 


.20209 


2451953 


.70209 


-0. 


09 


± 





.03 


2001 


Feb 


13. 


.21325 


2451953 


.71325 


-0. 


05 


± 





.03 


2001 


Feb 


13. 


.22439 


2451953 


.72439 


-0. 


15 


± 





.03 


2001 


Feb 


13. 


.23554 


2451953 


.73554 


-0. 


19 


± 





.0 1 


2001 


Feb 


13. 


.24671 


2451953 


.74671 


0. 


00 


± 





.01 


2001 


Feb 


14. 


.13972 


2451954 


.63972 


-0. 


05 


± 





.02 


2001 


Feb 


14. 


.15104 


2451954 


.65104 


-0. 


12 


± 





.02 


2001 


Feb 


14. 


.16228 


2451954 


.66228 


-0. 


25 


± 





.02 


2001 


Feb 


14. 


.17347 


2451954 


.67347 


-0. 


18 


± 





.02 


2001 


Feb 


14. 


.18477 


2451954 


.68477 


-0. 


11 


± 





.03 


2001 


Feb 


14. 


.19600 


2451954 


.69600 


-0. 


05 


± 





.03 


2001 


Feb 


14. 


.20725 


2451954 


.70725 


0. 


00 


± 





.03 


2001 


Feb 


14. 


.21860 


2451954 


.71860 


0. 


03 


± 





.03 


2001 


Feb 


14. 


.22987 


2451954 


.72987 


0. 


11 


± 





.04 


2001 


Feb 


14. 


.24112 


2451954 


.74112 


0. 


21 


± 





.04 


2001 


Feb 


14. 


.25234 


2451954 


.75234 


0. 


20 


± 





.05 


2001 


Feb 


14. 


.26356 


2451954 


.76356 


0. 


16 


± 





.05 


2001 


Feb 


15. 


.14518 


2451955 


.64518 


-0. 


06 


± 





.05 


2001 


Feb 


15. 


.15707 


2451955 


.65707 


-0. 


08 


± 





.02 


2001 


Feb 


15. 


.16831 


2451955 


.66831 


-0. 


15 


± 





.05 


2001 


Feb 


15. 


.19086 


2451955 


.69086 


0. 


05 


± 





.06 


2001 


Feb 


15. 


.20234 


2451955 


.70234 


0. 


19 


± 





.07 


2001 


Feb 


15. 


.23127 


2451955 


.73127 


0. 


01 


± 





.05 



a UT date at the start of the exposure 
^Julian date at the start of the exposure 
c Mean-subtracted apparent red magnitude; errors in- 
clude uncertainties in relative and absolute photometry 
added quadratically 



IPeixinho et"aTl pOOl hereafter |PDR 02) observed this 
object in September 2000, but having only one night's 
worth of data, they did not succeed in determining this 
object's rotational period unambiguously. We used their 
data to solve the degeneracy in our PDM result. The 
PDR02 data have not been absolutely calibrated, and 
the magnitudes are given relative to a bright field star. 
To be able to combine it with our own data we had 
to subtract the mean magnitudes. Our periodogram 
of (35671) 1998SNi65 (centered on the broad minimum) 
is shown in Fig. \S]p and can be compared with the re- 
vise d period ogram obtained with our data combined with 
the |PDR03 data (Fig. Et). The minima become much 
clearer with the additional data, but because of the 1- 
year time difference between the two observational cam- 
paigns, many close aliases appear in the periodogram. 
The absolute minimum, at P = 8.84 hr, corresponds to 
a double peaked lightcurve (see Fig. 0} . The second best 
fit, P — 8.7 hr, produces a more s cattered phase plot, 
in which the peak in the PDR02 data coincides with 
our night 2. Period P = 8.84 hr was also favored by 
the Monte Carlo method described in Section l3~2l being 
identified as the best fit in 55% of the cases versus 35% 
for P = 8.7 hr. The large size of the error bars compared 
to the amplitude is responsible for the ambiguity in the 
result. We will use P = 8.84 hr in the rest of the paper 
because it was consistently selected as the best fit. 

The amplitude, obtained using the Monte Carlo 
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TABLE 8 

Relative photometry measurements of 
2001 CZ31 



-0.2 



-0.1 









& 0.1 


UT Date a 


Julian Date b 


[mag] 


0.2 



2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 
2001 



Feb 28. 
Feb 28. 
Feb 28. 
Feb 28. 
Feb 28. 
Feb 28. 
Feb 28. 
Mar 01 
Mar 01 
Mar 01 
Mar 01 
Mar 01 
Mar 01 
Mar 01 
Mar 01 
Mar 03 
Mar 03 
Mar 03 
Mar 03 
Mar 03 
Mar 03 
Mar 03 



92789 
93900 
95013 
96120 
97235 
98349 
99475 
.00706 
.01817 
.02933 
.04046 
.05153 
.06304 
.08608 
.09808 
.01239 
.02455 
.03596 
.04731 
.05865 
.07060 
.08212 



2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451969. 
2451971. 
2451971. 
2451971. 
2451971. 
2451971. 
2451971. 
2451971. 



42789 
43900 
45013 
46120 
47235 
48349 
49475 
50706 
51817 
52933 
54046 
55153 
56304 
58608 
59808 
51239 
52455 
53596 
54731 
55865 
57060 
58212 



0.03 ± 0.05 
0.06 ± 0.04 
0.03 ± 0.04 
-0.09 ± 0.04 
-0.10 ± 0.04 
-0.12 ± 0.04 
-0.14 ± 0.03 



-0.02 
0.00 
0.03 
0.07 
0.10 
0.06 



0.03 
0.03 
0.03 
0.04 
0.04 
0.04 



± 
± 
± 
± 
± 
± 

-0.05 ± 0.04 
-0.05 ± 0.05 

0.15 ± 0.05 
-0.01 ± 0.05 

0.00 ± 0.04 
-0.02 ± 0.03 
-0.08 ± 0.04 
-0.04 ± 0.04 

0.01 ± 0.03 



a UT date at the start of the exposure 
Julian date at the start of the exposure 

c Mcan-subtractcd apparent red magnitude; errors in- 
clude uncertainties in relative and absolute photometry 
added quadratically 




10.0 8.5 
Period (ho 



Fig. 3.— Periodogram for the data on (35671) 1998SNi 65 . The 
lower left panel (b) shows an enlarged section near the minimum 
calculated using only the data published in this paper, and the 
lower right panel (c) shows the same region recalculated after 
adding the PDR02 data. 



0.2 



0.4 0.6 
Phase 



0.8 



Fig. 4.— Lightcurve of (35671) 1998SNi 65 . The figure represents 
the data phased with the best fit period P = 8.84 hr. Different 
symbols correspond to different nights of observation. The gray 
line is a 2nd order Fourier series fit to the data. The PDR02 data 
are shown as crosses. 




10 15 
Period (hr) 

Fig. 5.— Periodogram for the (79983) 1999 DF 9 data. The min- 
imum corresponds to a lightcurve period P = 6.65 hr. 




0.4 0.6 
Phase 

Fig. 6.— Same as Fig. E| for KBO (79983) 1999 DF9. The best 
fit period is P = 6.65 hr. The lines represent a 2nd order (solid 
line) and 5th (dashed line) order Fourier series fits to the data. 
Normalized \ 2 values of the fits are 2.8 and 1.3 respectively. 

method described in Section 13.31 is Am = 0.16 ± 
0.01 mag. This value was calculated using only our 
dat a, but i t did not change when recalculated adding 
the lPDR02l data. 

4.2. 1999 DF 9 

(79983) 1999 DFg shows large amplitude photometric 
variations (ArriR ~ 0.4 mag). The PDM periodogram 
for (79983) 1999 DFg is shown in Fig. The best-fit pe- 
riod is P = 6.65 hr, which corresponds to a double-peak 
lightcurve (Fig.[SJ). Other PDM minima are found close 
to P/2 « 3.3 hr and 9.2 hr, a 24 hr alias of the best pe- 
riod. Phasing the data with P/2 results in a worse fit 
because the two minima of the double peaked lightcurve 
exhibit significantly different morphologies (Fig. |5J; the 
peculiar sharp feature superimposed on the brighter min- 
imum, which is reproduced on two different nights, may 
be cau sed by a non-convex feature on the surface of the 
KBO llTorppaet al J 12003]) . Period P = 6.65 hr was se- 
lected in 65 of the 100 Monte Carlo replications of the 
dataset (see Section l3~2*|) . The second most selected so- 
lution (15%) was at P — 9hr. We will use P = 6.65 hr 
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TABLE 9 

Lightcurve Properties of Observed KBOs 



Object Designation 




Nts 


Amu c 


d 
1 




[mag] 




[mag] 


n 1 
[hr] 


(35671) 1998SNi B5 


21.20±0.05 


2(1) 


0.16±0.01 


8.84 (8.70) 


(79983) 1999 DF 9 




3 


0.40±0.02 


6.65 (9.00) 


2001 CZ31 




2(1) 


0.21±0.02 


4.71 (5.23) 


(19308) 1996 TOee 


21.26±0.06 


1 


? 




1996 TS 6 6 


21.76±0.05 


3 


<0.15 




(19521) Chaos 


20.74±0.06 


2 


<0.10 




(80806) 2OOOCM105 




2 


<0.14 




(66652) 1999RZ253 




3 


<0.05 




(47171) 1999 TC 36 




3 


<0.07 




(38628) Huya 




2 


<0.04 





a ncan red magnitude. Errors include uncertainties in relative and abso- 
lute photometry added quadratically; 

^number of nights with useful data. Numbers in brackets indicate num- 
ber of nights of data from other obser vers used for pe riod determination. 
Data for (35671) 1998 SNi 65 taken from lPeixinho et alJ 120021) and data for 
2001 CZ31 taken from SJ02; 

c lightcurvc amplitude; 

^lightcurve period (values in parenthesis indicate less likely solutions not 
entirely ruled out by the data). 




5.2 4.6 
Period (ho 

Fig. 7.— Periodogram for the 2001 CZ31 data. The lower left 
panel (b) shows an enlarged section near the minimum calculated 
using only the data published in this paper, and the lower right 
panel (c) shows the same region recalculated after adding the SJ02 
data. 



for the rest of the paper. 

The amplitude of the lightcurve, estimated as 
scribed in Section 1X51 is Am^ = 0-40 ± 0.02 mag. 



de- 



4.3. 2001 CZai 

This object shows substantial brightness variations 
(4.5cr above the comparison stars) in a systematic man- 
ner. The first night of data seems to sample nearly one 
complete rotational phase. As for (35671) 1998 SN165, 
the 2001 CZ31 data span only two nights of observations. 




0.4 0.6 
Phase 

Fig. 8.— Same as Fig. g] for KBO 2001 CZ31. The data are 
phased with period P = 4.71 hr. The points represented by crosses 
are taken from SJ02. 

In this case, however, the PDM minima (see Figs. [7^, 
and [3)) are very narrow and only two correspond to in- 
dependent periods, P = 4.69 hr (the minimum at 5.82 is 
a 24 hr alias of 4.69 hr), and P = 5.23 hr. 

2001 CZ31 has also been observed bv lSJ02l in February 
and April 2001, with inconclusive results. We used their 
data to try to rule out one (or bot h) of the two periods 
we found. We mean-subtracted the lSJ02l data in order to 
combine it with our uncalibrated observations. FigureEt 
shows the section of the periodogram around P = 5 hr, 
recalculated using ISJ02T s first night plus our own data. 
The aliases are due to the 1 month time difference be- 
tween the two observing runs. The new PDM minimum 
is at P — 4.71 hr - very close to the P = 4.69 hr deter- 
mined from our data alone. 

Visual inspection of the combined data set phased with 
P = 4.71 hr shows a very good match between lSJOffs firs t 
night (2001 Feb 20) and our own data (see Fig.®. SJ02's 
second and third nights show very large scatter and were 
not included in our analysis. Phasing the data with P = 
5.23 hr yields a more scattered lightcurve, which confirms 
the PDM result. The Monte Carlo test for uniqueness 
yielded P = 4.71 hr as the best-fit period in 57% of the 
cases, followed by P — 5.23 hr in 21%, and a few other 
solutions, all below 10%, between P = 5 hr and P = 6 hr. 
We will use P — 4.71 hr throughout the rest of the paper. 

We measured a lightcurve amp litude of Am = 0.21 ± 
0.02 mag. If we use both ours an d lSJ02f s first night data, 
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Am rises to 0.22 mag. 

4.4. Flat Lightcurves 

The fluctuations detected in the optical data on KBOs 
(19308) 1996 T0 66 , 1996 TS 66 , (47171) 1999 TC 36 , 
(66652) 1999 RZ253, (80806) 2000 CMi 05 , and 
(38628) Huya, are well within the uncertainties. 
(19521) Chaos shows some variations with respect to 
the comparison stars but no period was found to fit all 
the data. See Table [5] and Fig. [5] for a summary of the 
results. 

4.5. Other lightcurve measurements 

The KBO lightcurve database has increased con- 
siderably in the last few years , largely due to the 
ob servational campaig n of | S J02L with recent updates 
in ISheppard fc Jewitti l|2003|) and ISheppard fe JewitH 
l)2004[) . These authors have published observations 
and rotational data for a total of 30 KBOs (their 
SJ02 paper includes data for 3 other previously 
published lightcurves in the analysis). Other re- 
cently published KBO lightcurve s include those for 
(50000) Quaoar llOrtiz et all I2003D and the scattered 
KBO r2998!U999TDm URousselot et al .1 120031 . Of the 
10 KBO lightcurves presented in this paper, 6 are new 
to the database, bringing the total to 41. 

Table ITU1 lists the rotational data on the 41 KBOs that 
will be analyzed in the rest of the paper. 

5. ANALYSIS 

In this section we examine the lightcurve properties 
of KBOs and compare them with those of main-belt 
asteroids (MB As). The lightcurve data for these two 
families of objects cover different size ranges. MBAs, 
being closer to Earth, can be observed down to much 
smaller sizes than KBOs; in general it is very difficult 
to obtain good quality lightcurves for KBOs with diam- 
eters D < 50 km. Furthermore, some KBOs surpass the 
1000 km barrier whereas the largest asteroid, Ceres, does 
not reach 900 km. This will be taken into account in the 
analysis. 

The lightcurve data for asteroids were taken from 
the Harris Lightcurve Catalog 3 , Version 5, while 
the diameter data were obtained from the Lowell 
Observatory database of asteroids orbital elements 4 . 
The sizes of most KBOs were calculated from their 
absolute magnitude assuming an albedo of 0.04. 
The exceptions are (47171) 1999 TC 36 , (38 6 38) Huya, 
(28978) Ixion, (55636) 2002 TX 36 , (66 6 52) 1999 RZ 36 , 
(26308) 1998 SMi 65 , and (20000) Varuna for which the 
albed o has been shown to be inconsistent with the value 
0.04 l)Crundv et alJl2005l) . For example, in the case of 
(20000) Varuna simultaneous thermal and optical obser- 
vations have yielded a red geometric albedo of 0.0701q;q 3 " 
(jJewitt et aLlExnj) . 

5.1. Spin period statistics 

As Fig.llOlshows. the spin period distributions of KBOs 
and MBAs are significantly different. Because the sam- 
ple of KBOs of small size or large periods is poor, to 

3 http: //pdssbn. astro .umd. edu/ sbnhtml/asteroids/colors_lc 

4 f tp : //f tp . lowell . edu/pub/elgb/astorb . html 



avoid bias in our comparison we consider only KBOs 
and MBAs with diameter larger than 200 km and with 
periods below 20 hr. In this range the mean rotational 
periods of KBOs and MBAs are 9.23 hr and 6.48 hr, re- 
spectively, and the two are different with a 98.5% confi- 
dence according to Student's i-test. However, the differ- 
ent means do not rule that the underlying distributions 
are the same, and could simply mean that the two sets 
of data sample the same distribution differently. This 
is not the case, however, according to the Kolmogorov- 
Smirnov (K-S) test, which gives a probability that the 
periods of KBOs and MBAs are drawn from the same 
parent distribution of 0.7%. 

Although it is clear that KBOs spin slower than as- 
teroids, it is not clear why this is so. If collisions have 
contributed as significantly to th e angular momentum 
of KBOs as they did for MBAs ijFarinella et alJ 119821 
iCatullo et al.lll984j) . then the observed difference should 
be related to how these two families react to collision 
events. We will address the question of the collisional 
evolution of KBO spin rates in a future paper. 

5.2. Lightcurve amplitudes and the shapes of KBOs 

The cumulative distribution of KBO lightcurve ampli- 
tudes is shown in Fig. 1111 It rises very steeply in the 
low amplitude range (Am < 0.15 mag), and then be- 
comes shallower reaching large amplitudes. In quantita- 
tive terms, ~ 70% of the KBOs possess Am < 0.15 mag, 
while ~ 12% possess Am > 0.40 mag, with the maximum 
value being Am = 0.68 mag. [Note: Fig. ^2 does not 
include the KBO 2001 QG298 which has a lightcurve am- 
plitude Am = 1.14 ±0.04 mag, and would further extend 
the range of amplitudes. We do not include 2001 QG298 
in our analysis because it is tho ught to be a contact bi- 
nary llSheppard fc Jewitti [2004)] . Figure 1111 also com- 
pares the KBO distribution with that of MBAs. The 
distributions of the two populations are clearly distinct: 
there is a larger fraction of KBOs in the low amplitude 
range (Am < 0.15 mag) than in the case of MBAs, and 
the KBO distribution extends to larger values of Am. 

FigurclT2lshows the lightcurve amplitude of KBOs and 
MBAs plotted against size. KBOs with diameters larger 
than D = 400 km seem to have lower lightcurve am- 
plitudes than KBOs with diameters smaller than D = 
400 km. Student's t-test confirms that the mean ampli- 
tudes in each of these two size ranges arc different at the 
98.5% confidence level. For MBAs the transition is less 
sharp and seems to occur at a smaller size (D ~ 200 km) . 
In the case of asteroids, the accepted explanation is 
that small bodies (D < 100 km) are fragments of high- 
velocity impacts, whereas of thei r larger counterpart s 
(D > 200 km) generally are not <|Catullo et, al] li984L 
The lightcurve data on small KBOs are still too sparse 
to permit a similar analysis. In order to reduce the effects 
of bias related to body size, we can consider only those 
KBOs and MBAs with diameters larger than 200 km. In 
this size range, 25 of 37 KBOs (69%) and 10 of 27 MBAs 
(37%) have lightcurve amplitudes below 0.15 mag. We 
used the Fisher exact test to calculate the probability 
that such a contingency table would arise if the lightcurve 
amplitude distributions of KBOs and MBAs were the 
same: the resulting probability is 0.8%. 
html The distribution of lightcurve amplitudes can be used 
to infer the shapes of KBOs, if certain reasonable as- 
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Fig. 9. — The "flat" lightcurves are shown. The respective amplitudes are within the photometric uncertainties. 



sumptions are made (see, e.g., lLL03aj) . Generally, ob- 
jects with elongated shapes produce large brightness vari- 
ations due to their changing projected cross-section as 
they rotate. Conversely, round objects, or those with the 
spin axis aligned with the line of sight, produce little or 
no brightness variations, resulting in "flat" lightcurves. 
Figure IP21 shows that the lightcurve amplitudes of KBOs 
with diameters smaller and larger than D = 400 km 



are significantly different. Docs this mean that the 
shapes of KBOs are also different in these two size 
ranges? To investigate this possibility of a size depen- 
dence among KBO shapes we will consider KBOs with 
diameter smaller and larger than 400 km separately. We 
shall loosely refer to objects with diameter D > 400 km 
and D < 400 km as larger and smaller KBOs, respec- 
tively. 
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TABLE 10 

Database of KBOs Lightcurve Properties 



Object Designation Class a Size b P c Am R d Ref. 

[km] [hr] [mag] 



KBOs considered to have Am < 0.15 mag 



(15789) 


1993 SC 


P 


240 


0.04 


7, 2 


(15820) 


1994 TB 


P 


220 


<0.04 


10 


(26181) 


1996 GQ 91 


S 


730 


<0.10 


10 


(15874) 


1996 TL 66 


s 


480 


0.06 


7, 4 


(15875) 


1996 TP 6 6 


p 


250 


0.12 


7, 1 


(79360) 


1997 CS29 


C 


630 


<0.08 


10 


(91133) 


1998 HK151 


p 


170 


<0.15 


10 


(33340) 


1998 VG44 


p 


280 


<0.10 


10 


(19521) 


Chaos 


C 


600 


<0.10 


13, 10 


(26375) 


1999 DEg 


s 


700 


<0.10 


10 


(47171) 


1999 TC36 


Pb 


300 


<0.07 


13, 11 


(38628) 


Huya 


P 


550 


<0.04 


13, 11 


(82075) 


2OOOYW134 


s 


790 


<0.10 


11 


(82158) 


2001FPis5 


s 


400 


<0.10 


11 


(82155) 


2OOIFZ173 


s 


430 


<0.06 


10 




2001 KD77 


p 


430 


<0.10 


11 


(28978) 


Ixion 


p 


820 


<0.10 


11, 5 




2001 QF 298 


p 


580 


<0.10 


11 


(42301) 


2001 UR163 


s 


1020 


<0.10 


11 


(42355) 


2002 CR 46 


s 


210 


<0.10 


11 


(55636) 


2002 TX300 


C 


710 16.24 


0.08±0.02 


11 


(55637) 


2002 UX 25 


C 


1090 


<0.10 


11 


(55638) 


2002 VE95 


p 


500 


<0.10 


11 


(80806) 


2OOOCM105 


C 


330 


<0.14 


1.3 


(66652) 


1999 RZ253 


Cb 


170 


<0.05 


13 




1996 TS 66 


C 


300 


<0.14 


13 



KBOs considered to have Am > 0.15 mag 



(32929) 


1995 QY 9 


P 


180 


7.3 


0.60±0.04 


10, 7 


(24835) 


1995 SM 55 


C 


630 


8.08 


0.19±0.05 


11 


(1930S) 


1996 TOee 


C 


720 


7.9 


0.26±0.03 


11, 3 


(26308) 


1998SMi 65 


R 


240 


7.1 


0.45±0.03 


10, 8 


(33128) 


1998 BU 4 8 


S 


210 


9.8 


0.68±0.04 


10, 8 


(40314) 


1999 KRi 6 


C 


400 


11.858 


0.18±0.04 


10 


(47932) 


2000 GN171 


p 


360 


8.329 


0.61±0.03 


10 


(20000) 


Varuna 


C 


980 


6.34 


0.42±0.03 


10 


2003 AZ 84 


p 


900 


13.44 


0.14±0.03 


11 




2001 QG 29g 


Pcb 


240 


13.7744 


1.14±0.04 


12 


(50000) 


Quaoar 


C 


1300 


17.6788 


0.17±0.02 


6 


(29981) 


1999 TD10 


s 


100 


15.3833 


0.53±0.03 


9 


(35671) 


1998SNi65 


C 


400 


8.84 


0.16±0.01 


13 


(79983) 


1999 DFg 


C 


340 


6.65 


0.40±0.02 


13 




2001 CZ31 


C 


440 


4.71 


0.21±0.06 


13 



Ref erences. (1) ICollander-Brown et all 1199911: (2 ) IDavies et"aT1 

119971): (3 ) IHainaut et al] 1200011: (4) ILuu fc Jewitti 1199811: (5) 

IQrtiz et all 1200111: (61 IQrtiz et al.l 1200311 * ( 7) IRomanishin fe Teglerl 
119991): Ts) IRom anishin e t alJ 1200111: (9) IR.ousselot et al.l 120031 : 
flO) IShep pard fc Jew*tJ T2*oTj*li : fill ISheppard fe Jewitti 12003ll : (12) 

ISheppard £^ewitTl2(*r*"lll : (13) this work. 

a Dynamical class (C — classical KBO, P — Plutino, R. — 2:1 Resonant b — 
binary KBO); 

Diameter in km assuming an albedo of 0.04 except when measured (see 
text); 

c Pcriod of the lightcurve in hours. For KBOs with both single and double 
peaked possible lightcurves the double peaked period is listed; 
^Lightcurve amplitude in magnitudes. 



We approximate the shapes of KBOs by triaxial ellip- 
soids with semi-axes a > b > c. For simplicity we con- 
sider the case where b — c and use the axis ratio a = a/b 
to characterize the shape of an object. The orientation 
of the spin axis is parameterized by the aspect angle 8, 
defined as the smallest angular distance between the line 
of sight and the spin vector. On this basis the lightcurve 
amplitude Am is related to a and 8 via the relation (Eq. 



(2) of [TJ,03al with c = 1) 

Am = 2.5 log - 



2 5 2 



v i + a 2 + (a 2 -i) cos(26>) ' W 

Following LL03b we model the shape distribution by a 
power-law of the form 



/(a) da oc a q da 



(2) 
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Fig. 10. — Histograms of the spin periods of KBOs (upper 
panel) and main belt asteroids (lower panel) satisfying D > 200 km, 
Am > 0.15 mag, P < 20 hr. The mean rotational periods of KBOs 
and MBAs are 9.23 hr and 6.48 hr, respectively. The y-axis in both 
cases indicates the number of objects in each range of spin periods. 
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Fig. 11. — Cumulative distribution of lightcurve amplitude for 
KBOs (circles) and asteroids (crosses) larger than 200 km. We 
plot only KBOs for which a period has been determined. KBO 
2OOIQG2981 thought to be a contact binary CShcpoard fc Jewitul 
120041) . is not plotted although it may be considered an extreme 
case of elongation. 
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Fig. 12. — Lightcurve amplitudes of KBOs (black circles) and 
main belt asteroids (gray crosses) plotted against object size. 

where f(a) da represents the fraction of objects with 
shapes between 5, and a + da. We use the measured 
lightcurve amplitudes to estimate the v alue of q by em- 
ploying both the method described in ILL03al and by 
Monte Carlo f itting the observed amplitude distribution 
l|SJ02l lLLb 3b). The latter consists of generating artifi- 
cial distributions of Am (Eq. ^) with values of a drawn 
from distributions characterized by different g's (Eq. |2J), 
and selecting the one that best fits the observed cumu- 
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Fig. 13. — Observed cumulative lightcurve amplitude distribu- 
tions of KBOs (black circles) with diameter smaller than 400 km 
(upper left panel), larger than 400 km (lower left panel), and of all 
the sample (lower right panel) are shown as black circles. Error 
bars are Poissonian. The best fit power-law shape distributions of 
the form f(a) da = a~ q da were converted to amplitude distribu- 



l s), a 

m 



shown as solid lines. The best fit q's are listed in Table [ 

lative amplitude distribution (Fig. Ill|) . The values of 6 
are generated assuming random spin axis orientations. 
We use the K-S test to compare the different fits. The 
errors ar e derived by bootstrap resampling the original 
data set (Efron1 ll979j) . and measuring the dispersion in 
the distribution of best-fit power-law indexes, qi, found 
for each bootstr ap repli cation. 

Following the ILL03al method we calculate the proba- 
bility of finding a KBO with Am > 0.15 mag: 



p(Am > 0.15) 



/(a). 



K 



(a 2 - l)K 



da. (3) 



where K = io°- 8 >< 015 , /(g) = Ca^, and C is 
a normalization constant. This probability is calcu- 
lated for a range of q's to determine the one that 
best matches the observed fraction of lightcurves with 
amplitude larger than 0.15 mag. These fractions are 
/(Am > 0.15mag;L> < 400km) = 8/19, and 
/(Am > 0.15 mag; D > 400 km) = 5/21, and /(Am > 
0.15 mag) = 13/40 for the complete set of data. The re- 
sults are summarized in Table ITT1 and shown in Fig. 1131 
The uncertainties in the values of q obtained using the 
lLL03al method (q = 4.3±^ for KBOs with D < 400 km 

and q = 7A±l^ for KBOs with D > 400 km ; see Ta- 
ble 111(1 do not rule out similar shape distributions for 
smaller and larger KBOs. This is not the case for the 
Monte Carlo method. The reason for this is that the 
LL03a method relies on a single, more robust parame- 
ter: the fraction of lightcurves with detectable variations. 
The sizeable error bar is indicative that a larger dataset is 
needed to better constrain the values of q. In any case, 
it is reassuring that both methods yield steeper shape 
distributions for larger KBOs, implying more spherical 
shapes in this size range. A distribution with q ~ 8 
predicts that ~75% of the large KBOs have a/6 < 1.2. 
For the smaller objects we find a shallower distribution, 
q ~ 4, which implies a significant fraction of very elon- 
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TABLE 11 
Best fit parameter to the KBO shape 

DISTRIBUTION 



Size Range '' 



Method b 
LL03 MC 



D < 400 km 
D > 400 km 
All sizes 



4 n+2.o 
*-°-1.6 
. 7 4 +3.l 

+ 1.6 
1.3 



: 5.7 



q = 3.8 ±0.8 
q = 8.0 ± 1.4 
g = 5.3 ±0.8 



Range of KBO diameters, in km, considered 
i each case; 

3 LL03 is the method described in 



ILacerda fc Luul 1200311 . and MC is a Monte 
Carlo fit of the lightcurve amplitude distribu- 
tion. 

gated objects: ~20% have a/b > 1.7. Although based on 
small numbers, the shape distribution of large KBOs is 
well fit by a simple power-law (the K-S rejection proba- 
bility is 0.6%). This is not the case for smaller KBOs for 
which the fit is poorer (K-S rejection probability is 20%, 
see Fig. 11 31) . Our results are in agreement with previous 
studies of the overall KBO shape distribution, which had 
already shown that a simple po wer-law does n ot explain 
the shapes of KBOs as a whole (LL03b, SJ02). 

The results presented in this section suggest that the 
shape distributions of smaller and larger KBOs are dif- 
ferent. However, the existing number of lightcurves is 
not enough to make this difference statistically signifi- 
cant. When compared to asteroids, KBOs show a pre- 
ponderance of low amplitude lightcurves, possibly a con- 
sequence of their possessing a larger fraction of nearly 
spherical objects. It should be noted that most of our 
analysis assumes that the lightcurve sample used is ho- 
mogeneous and unbiased; this is probably not true. Dif- 
ferent observing conditions, instrumentation, and data 
analysis methods introduce systematic uncertainties in 
the dataset. However, the most likely source of bias in 
the sample is that some flat lightcurves may not have 
been published. If this is the case, our conclusion that 
the amplitude distributions of KBOs and MBAs are dif- 
ferent would be strengthened. On the other hand, if most 
unreported non-detections correspond to smaller KBOs 
then the inferred contrast in the shape distributions of 
different-sized KBOs would be less significant. Clearly, 
better observational contraints, particularly of smaller 
KBOs, are necessary to constrain the KBO shape distri- 
bution and understand its origin. 

5.3. The inner structure of KBOs 

In this section we wish to investigate if the rotational 
properties of KBOs show any evidence that they have 
a rubble pile structure; a possible dependence on ob- 
ject size is also investigated. As in the case of aster- 
oids, collisional evolution may have played an important 
role in modifying the inner structure of KBOs. Large 
asteroids (D > 200 km) have in principle survived col- 
lisional destruction for the age of the solar system, but 
may nonetheless have been converted to rubble piles by 
repeated impacts. As a result of multiple collisions, the 
"loose" pieces of the larger asteroids may have reassem- 
bled into shapes clos e to triaxial equilibrium ellipsoids 
ijFarinella et al.l[l981j) . Instead, the shapes of smaller 
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Fig. 14. — Lightcurve amplitudes versus spin periods of KBOs. 
The black filled and open circles represent objects larger and 
smaller than 400 km, respectively. The smaller gray circles show 
the results of numeri cal simulations of "rubble-pile" collisions 
(Lcinhardt ct al. 2000). The lines represent the locus of rotat- 
ing ellipsoidal figures of hydrostatic equilibrium with densities 
p = 500 kg m -3 (dotted line), p = 1000 kg m — 3 (shorter dashes) 
and p = 2000 kg m~ 3 (longer dashes). Shown in grey next to the 
line s are the axis ratios, a /b, of the ellipsoidal solutions. Both 
the ILeinhardt e~ al. (2000) results and the figures of equilibrium 
assume equator-on observing geometry and therefore represent up- 
per limits. 

asteroids ( D < 100 km) are c onsistent with collisional 
fragments ijCatullo et al.lll98^1 . indicating that they are 
most likely by-products of disruptive collisions. 

Figure 1141 plots the lightcurve amplitudes versus spin 
periods for the 15 KBOs whose lightcurve amplitudes 
and spin period are known. Open and filled symbols in- 
dicate the KBOs with diameter smaller and larger than 
D = 400 km, respectively. Clearly, the smaller and 
larger KBOs occupy different regions of the diagram. 
For the larger KBOs (black filled circles) the (small) 
lightcurve amplitudes are almost independent of the ob- 
jects' spin periods. In contrast, smaller KBOs span a 
much broader range of lightcurve amplitudes. Two ob- 
jects have very low amplitudes: (35671) 1998SNi 65 and 
1999 KRi6, which have diameters D ~ 400 km and fall 
precisely on the boundary of the two size ranges. The 
remaining objects hint at a trend of increasing Am with 
lower spin rates. The one exception is 1999 TDio, a 
Scattered Disk Object (e = 0.872, a = 95.703 AU) that 
spends most of its orbit in rather empty regions of space 
and most likely has a different collisional history. 

For comparison, Fig. ^] also shows results of N- 
body simulations of collisi ons between "idea l" rubble 
piles (gray filled circles: ILeinhardt et al.ll2000[K and the 
lightcurve amplitude-spin period relation predicted by el- 
lipsoidal fig ures of hydrostatic equilibrium ( dashed and 
dotted lines: lChandrasekharlll969l lHolsappldl200l . The 
latter is calculated from the equilibrium shapes that ro- 
tating uniform fluid bodies assume by balancing grav- 
itational and centrifugal acceleration. The spin rate- 
shape relation in the case of uniform fluids depends solely 
on the density of the body. Although fluid bodies be- 
have in many respects differently from rubble piles, they 
may, as an extreme case, provide insight on the equi- 
librium shapes of gravitationally bound agglomerates. 
The lightcurve amplitudes of both theoretical expecta- 
tions are plotted assuming an equator-on observing ge- 
ometry. They should therefore be taken as upper limits 
when compared to the observed KBO amplitudes, the 
lower limit being zero amplitude. 

The simulations of ILeinhardt et alJ <|2000t hereafter 



14 



Lacerda & Luu 



LRQ) consist of collisions between agglomerates of small 
spheres meant to simulate collisions between rubble piles. 
Each agglomerate consists of ~ 1000 spheres, held to- 
gether by their mutual gravity, and has no initial spin. 
The spheres are indestructible, have no sliding friction, 
and have coefficients of restitution of ~ 0.8. The bulk 
density of the agglomerates is 2000 kg m -3 . The impact 
velocities range from ~ zero at infinity to a few times the 
critical velocity for which the impact energy would ex- 
ceed the mutual gravitational binding energy of the two 
rubble piles. The impact geometries range from head-on 
collisions to grazing impacts. The mass, final spin period, 
and shape of the largest remn ant of each collision are 
registered (see Table 1 of LRQ). From their results, we 
selected the outcomes for which the mass of the largest 
remnant is equal to or larger than the mass of one of the 
colliding rubble piles, i.e., we selected only accreting out- 
comes. The spin periods and lightcurve amplitudes that 
would be generated by such remnants (assuming they 
are observed equator-on) are plotted in Fig. ^] as gray 
circles. Note that, although the simulated rubble piles 
have radii of 1 km, since the effects of the collision scale 
with the ratio of impact energy to gravitational bindin g 
energy of the colliding bodies ijBenz fc Asphaud Fl999'). 
the mode l results should apply to other sizes. Clearly, 
the LRQ model makes several specific assumptions, and 
represents one possible idealization of what is usually re- 
ferred to as "rubble pile". Nevertheless, the results are 
illustrative of how collisions may affect this type of struc- 
ture, and are useful for comparison with the K BO da ta. 

The lightcurve amplitudes resulting from the LRQ| ex- 
periment are relatively small (Am < 0.25 mag) for spin 
periods larger than P ~ 5.5 hr (see Fig. I14fl . Objects 
spinning faster than P = 5.5 hr have more elongated 
shapes, resulting in larger lightcurve amplitudes, up to 
0.65 magnitudes. The latter are the result of collisions 
with higher an gular m omentum transfer than the former 
(see Table 1 of LRQ]) . The maximum spin rate attained 
by the rubble piles, as a result of the collision, is ~ 4.5 hr. 
This is consistent with the maximum spin expected for 
bodies in hydrostatic equilibrium with the same density 
as the rubble piles (p = 2000 k g m~ 3 - see long-dashed 
line in Fig. 113)1 . The results of LRQ: show that colli- 
sions between ideal rubble piles can produce elongated 
remnants (when the projectile brings significant angular 
momentum into the target), and that the spin rates of 
the collisional remnants do not extend much beyond the 
maximum spin permitted to fluid uniform bodies with 
the same bulk density. 

The distribution of KBOs in Fig. El is less clear. In- 
direct estimates of KBO bulk densities indicate values 
p ~ 1000 kg m" 3 lILuu fc JewitU l200l . If KBOs are 
strengthless rubble piles with such low densities we would 
not expect to find objects with spin periods lower than 
P ~ 6hr (dashed line in Fig. 114(1 . However, one ob- 
ject (2001 CZ31) is found to have a spin period below 
5hr. If this object has a rubble pile structure, its den- 
sity must be at least ~ 2000 kg m -3 . The remaining 
14 objects have spin periods below the expected upper 
limit, given their estimated density. Of the 14, 4 ob- 
jects lie close to the line corresponding to equilibrium 
ellipsoids of density p = 1000 kg m -3 . One of these 
objects, (20000) Va r una, has been studied in detail by 
ISheppard fc Jewittl (|2002fl . The authors conclude that 



(20000) Varuna is best interpreted as a rotationally de- 
formed rubble pile with p < 1000 kg m -3 . One object, 
2001 QG298, has an exceptionally large lightcurve ampli- 
tude (Am = 1.14 mag), indicative of a very elongated 
shape (axes ratio a/b > 2.85), but given its modest spin 
rate (P = 13.8 hr) and approximate size (D ~ 240 km) 
it is unlikely that it would be able to keep such an elon- 
gated shape against the c rush of gravity. Analysis o f the 
lightcurve of this object ((Sheppard fc JewiT3l2004() sug- 
gests it is a close/contact binary KBO. The same applies 
to two other KBOs, 2000 GN m and (33128) 1998 BU 48 , 
also very likely to be contact binaries. 

To summarize, it is not clear that KBOs have a rubble 
pile structure, based on their available rotational proper- 
ties. A comparison with computer simulations of rubble 
pile collisions shows that larger KBOs (D > 400 km) 
occupy the s ame r egion of the period-amplitude dia- 
gram as the ILRQI results. This is not the case for 
most of the smaller KBOs (D < 400 km), which tend 
to have larger lightcurve amplitudes for similar spin pe- 
riods. If most KBOs are rubble piles then their spin 
rates set a lower limit to their bulk density: one ob- 
ject (2001 CZ31) spins fast enough that its density must 
be at least p ~ 2000 kg m -3 , while 4 other KBOs (in- 
cluding (20000) Varuna) must have densities larger than 
p ~ 1000 kg m -3 . A better assessment of the inner struc- 
ture of KBOs will require more observations, and detailed 
modelling of the collisional evolution of rubble-piles. 

6. CONCLUSIONS 

We have collected and analyzed R-band photometric 
data for 10 Kuiper Belt objects, 5 of which have not 
been studied before. No significant brightness varia- 
tions were detected from KBOs (80806) 2000 CMi 05) 
(66652) 1999 RZ 253 , 1996 TS 66 . Previously ob- 

served KBOs (19521) Chaos, (47171) 1999 TC 36 , and 
(38628) Huya were confirmed to have very low ampli- 
tude lightcurves (Am < 0.1 mag). (35671) 1998 SNi 65 , 
(79983) 1999 DF 9 , and 2001 CZ31 were shown to have 
periodic brightness variations. Our lightcurve amplitude 
statistics are thus: 3 out of 10 (30%) observed KBOs 
have Am > 0.15 mag, and 1 out of 10 (10%) has 
Am > 0.40 mag. This is consistent with previously 
published results. 

The rotational properties that we obtained were com- 
bined with existing data in the literature and the total 
data set was used to investigate the distribution of spin 
period and shapes of KBOs. Our conclusions can be 
summarized as follows: 

1. KBOs with diameters D > 200 km have a mean 
spin period of 9.23 hr, and thus rotate slower on 
average than main belt asteroids of similar size 
((-P)mbAs = 6.48 hr). The probability that the two 
distributions are drawn from the same parent dis- 
tribution is 0.7%, as judged by the KS test. 

2. 26 of 37 KBOs (70%, D > 200 km) have lightcurve 
amplitudes below 0.15 mag. In the asteroid belt 
only 10 of the 27 (37%) asteroids in the same size 
range have such low amplitude lightcurves. This 
difference is significant at the 99.2% level according 
to the Fisher exact test. 

3. KBOs with diameters D > 400 km have lightcurves 
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with significantly (98.5% confidence) smaller am- 
plitudes ((Am) = 0.13 mag, D > 400km) than 
KBOs with diameters D < 400km ((Am) = 
0.25 mag, D < 400 km). 

4. These two size ranges seem to have different shape 
distributions, but the few existing data do not ren- 
der the difference statistically significant. Even 
though the shape distributions in the two size 
ranges are not inconsistent, the best-fit power- 
law solutions predict a larger fraction of round 
objects in the D > 400 km size range (f(a/b < 
1.2) - 70t\l%) than in the group of smaller ob- 
jects (f(a/b < 1.2) ~ 42±?g%). 

5. The current KBO lightcurve data are too sparse to 
allow a conclusive assessment of the inner structure 
of KBOs. 

6. KBO 2001 CZ 3 i has a spin period of P = 4.71 hr. 
If this object has a rubble pile structure then its 
density must be p > 2000 kg m -3 . If the object has 



a lower density then it must have internal strength. 

The analysis presented in this paper rests on the as- 
sumption that the available sample of KBO rotational 
properties is homogeneous. However, in all likelihood 
the database is biased. The most likely bias in the sam- 
ple comes from unpublished flat lightcurves. If a signifi- 
cant fraction of flat lightcurves remains unreported then 
points 1 and 2 above could be strengthened, depending 
on the cause of the lack of brightness variation (slow 
spin or round shape). On the other hand, points 3 and 4 
could be weakened if most unreported cases correspond 
to smaller KBOs. Better interpretation of the rotational 
properties of KBOs will greatly benefit from a larger and 
more homogeneous dataset. 
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